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ABSTRACT 

Context. High-precision frequencies of acoustic modes in red giant stars are now available thanks to the long observing length and 
high-quality of the light curves provided by the NASA Kepler mission, thus allowing to probe the interior of evolved cool low-mass 
stars with unprecedented level of detail. 

Aims. We characterize the acoustic signature of the helium second ionization zone in a sample of 18 low-mass low-luminosity red 
giants by exploiting new mode frequency measurements derived from more than four years of Kepler observations. 

Methods. We analyze the second frequency differences of radial acoustic modes in all the stars of the sample by using the Bayesian 
code Diamonds. 

Results. We find clear acoustic glitches due to the signature of helium second ionization in all the stars of the sample. We measure 
the acoustic depth and the characteristic width of the acoustic glitches with a precision level on average around ~2 % and ~8 %, 
respectively. We find good agreement with theoretical predictions and existing measurements from the literature. Lastly, we derive the 
amplitude of the glitch signal at v^ax for the second differences and for the frequencies with an average precision of ~6 %, obtaining 
values in the range 0.14-0.24//Hz, and 0.08-0.33//Hz, respectively, which can be used to investigate the helium abundance in the 
stars. 
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1. Introduction 


The so-called acoustic glitches are regions of sharp-structure 
variation located in the interior of the stars and caused by the 
presence of a change in the energy transport from radiative to 
convective, by a rapid variation in the chemical composition, or 
by ionization zones of chemical elements such as hydrogen and 
helium. As or iginally predicted for the Sun (e.g. Vorontsov 1988[ 
|Gough|1990| ), these regions produce tiny and regular variations 
in the frequency of the acoustic (p) modes that can be detected 
by direct measurement of the characteristic large frequency sep¬ 
aration, namely the frequency separation between modes having 
the same angular degree. 

By studying the glitch signature in the Sun, it was possible 
to measure the acoustic position of the base of the convective 
zone and of the helium second ionization (Heii) zone, as well 
as to provide estimates of the helium abundance in the envelope, 
and the extent of the overshooting (e.g.|Basu & Antia|1995"}|Bas u| 
1997t [Monteiro & Thompson|2QQ5[ |Cliriste nsen -Dalsgaard et al. 


2011). Expected also for d istant stars (e.g.|Monteiro et al.|20Q0[ 


Mazumdar & Antia|2001t [Ballot et al.|2004| ), thanks to the ad¬ 

vent o f the CoRoT ([Baglin et al.||2006|) and Kepler space mis¬ 
sions ( |Borucki et al.pO^lQt |Koch et aLpOlO} , which have re¬ 
leased an outstanding amount of high-quality photometric obser- 
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vations for thousands of stars, these frequency shifts have been 
discovered and analyzed in ma ny low-mass main-sequence, sub 
giant and red giant stars (RGs) (|Mosser et al.|2()Tm|Miglio et al 


20 10[ [Mazumdar et al.|2012[|2Q14[|Verma et al.|2Qi4| ), allowing 

to constrain both the position of the base of the convective zone 
and that of the He ii zone in main-sequence and sub-giant stars, 
and the position of the He ii zone in the case of the red giants. 

The asteroseismology of red giant stars, in particular, has 
brought several important breakthroughs in the stellar physics of 


low-mass stars in the latest years (e.g. Beck et al.|2011t Mosser 

et al. 

2011al|Bedding et al. 

20111 Beck et al.||2012| Deheuvels 


et al.|2012] r The characterization of the glitch signatures is able 

to provide tighter constraints on the chemical composition and 
the internal stratification of the star, and potentially allows to 
retrieve helium abunda nces in distant stars, ess ential for pop¬ 
ulation studies (e.g. see Broomhall et al. 2014| hereafter B14, 
and references therein). More recent studies focusing toward 
these evolved cool stars have analyzed the glitches due to the 
He II zone for an ensemble of more than a hundred targets ob¬ 
served by Kepler ( jVrard et al.||2014 ), and thoroughly investi¬ 
gated the properties of the signature from a theoretical point of 
view (B14, see also Christensen-Dalsgaard et al.|2014 for more 
discussion). 


The recent availability of Kepler datasets spanning more than 
four years of nearly continuous observations, coupled with the 
development of new computational advances in asteroseismic 
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data a nalysis (e.g. see Corsaro & De Ridder|2014 Corsaro et al.| 
2015[ hereafter C15), enables the study of the acoustic glitch 


signatures in red giants with an unprecedented level of detail. 

In this paper, we report on the evidence of clear acoustic 
glitch signatures due to the Hen zone in the sample of 19 red 
giants recently investigated by Cl5, hence we fully characterize 
the oscillatory signal by means of a Bayesian approach. 


2. Data analysis 

As noticed by Cl5, the low-mass low-luminosity red giants 
(LRGs) are well suited candidates to test stellar structure mod¬ 
els and stellar evolution theory. The less-evolved stage in the red 
giant branch (RGB) of the stellar evolution for the LRGs, im¬ 
plies the highest frequency of maximum power v^ax for a red 
giant (between 100 and 200//Hz), hence a broader power ex¬ 
cess caused by the oscillations, and consequently a larger num¬ 
ber of radial orders observed (in general between six and nine). 
By having a larger number of high signal-to-noise ratio p-mode 
frequencies available, one is thus able to constrain the signature 
of the glitches more efficiently. 

In this work we analyze the sample of LRGs studied by Cl5, 
who peak bagged their full oscillation spectrum using Kepler 
observations from QO till Q17.1, a total of -1470 days, with a 
frequency resolution of ^Vbm - 0.008 yuHz. The stars have Vmax 
values ranging from 110 to - 190yuHz and masses in the inter¬ 
val 1-2 Mo. We adopt the high-precision individual frequency 
measurements from C15, and follow the theoretical approach by 
B14. 

In the present analysis we refer to the first (frequency) differ¬ 
ence, as the large frequency separation of a given angular degree, 
as a function of the frequency in the power spectral density of the 
star, Av^ (v). For a radial order n, Av^ (v) is thus defined as 


= yn+U-yn,£- ( 1 ) 

where Vn,e is the central frequency of the mode with angular de¬ 
gree i and radial order n. In a ddition, we co mpute the second 
(frequency) difference (see e.g. |Gough||1990| ), A 2 y^ (v), defined 
for a single radial order as 

A2V„/ = Vn+U - + Vn-U 

= - Ay„_i,^ . 


W e fit the acoustic glitch sig natures with the model introduced 
by |Houdek & Gough] ( |2007| ), and used by B14 for RGs, defined 


as 


^2^0n4 = Aa)n,{ cos [2 (ths + ^)]+c, (3) 

with = 271 V and lS.20Jn,£ = 2;7rA2y„/, A a dimensionless 
amplitude of the signature signal, thc the acoustic depth of the 
He II zone, b its characteristic width, 0 and c a constant phase 
shift and offset, respectively, of the oscillatory signal. Following 
the arguments discussed by B14, we apply the fit to the second 
differences only because they are less prone to additional vary¬ 
ing components such as hydrogen ionization and non-adiabatic 
processes, and to the general frequency dependence of the large 
separation caused by the development of the se cond-order term 
of the asymptotic relation ( [Mosser et al.|2Ql lb] ). The second dif¬ 
ferences are at the same time available in a reasonably high num¬ 
ber of measurements (two less than the total number of modes 
obtained for a given angular degree), still allowing to constrain 
the model parameters without leading to degeneracies in the so¬ 
lutions. 


Despite the possible presence of the oscillatory component 
in modes of angular degree ^ > 0, we point out that only ra¬ 
dial mode frequencies are used for the final fit. The reason be¬ 
hind this choice mainly is the need to exploit pure p-mode char¬ 
acter oscillations (see also B14), which in the case of RGs are 
only available in the form of ^ = 0 modes. This is because the 
coupling occurring between p modes of angular degree ^ > 0 
and g modes arising from the radiative interior can hamper the 
asymptotic behavior of the corresponding modes by producing 
so-called mixed modes, whose frequencies deviate from the ex- 
pected position of a pure p mode oscillation (e.g. |Beck et al 

[20TT] ). 

We perform all the fi ts following a Bayesian app roach done 
by means of Diamonds ( [Corsaro & De Ridder|2014| ), hence ex¬ 
ploiting a nested sampling Monte Carlo method to perform the 
inference and estimate the free parameters of the model given 
by Eq. ^ from their individual marginal probability distribu¬ 
tions (see [Corsaro & De Ridder||2014| for more details). The 
configuring parameters of Diamonds (following the definitions 
by [Corsaro & De Ridder[[2014[ ) used for all the computations 
are: initial enlargement fraction 1.0 < /o < 1.7, shrinking rate 
a = 0.02-0.03, number of live points Nh^q = 1000, number 
of clusters 1 < A^ciust ^ 4, number of total drawing attempts 
4^attempts = 10"^, number of nested iterations before the first clus¬ 
tering Minit = 1000, and number of nested iterations with the 
same clustering Msame = 50. 

For this analysis we adopt a normal likelihood function, as 
that used by [Corsaro et al.[p013[ ), which takes into account the 
uncertainties, with correlations included, on the measurements 
of the second differences. This assumes that the residuals arising 
from the difference between predicted and measured second dif¬ 
ferences are Gaussian distributed. Concerning the set up of pri¬ 
ors, since we do not have initial guesses available from the litera¬ 
ture for the given stars, we use uniform (i.e. fiat) prior probability 
distributions for all the free parameters of Eq. ([^, with lower and 
upper boundaries for each parameter range obtained by compar- 
is on with existing measurements of the acoust ic depths derived 
by [Miglio et al. ( [2010[ ), [Mazumdar et al.[ ( [2014] ), and the theoret¬ 
ical results by B14 in the observed range of Vmax- The choice of 
uniform priors also yields a faster computation with Diamonds, 
as already discussed by Corsaro & De Ridder ( 2014[) and C15. 

Following the discussion by [Ballot et al.[ ( 2004[ ) and B14, 
we compute the acoustic radius of the Hen zones, ffieii, since 
it represents a quasi unbiased estimator of the acoustic posi¬ 
tion of the glitch. This is done by using the mean large fre¬ 
quency separation, (Av), obtained from the radial mode frequen¬ 
cies provided by C15, giving the total acoustic radius of the 
star, T = (2(Ay))“\ hence the acoustic radius of the Hen zone, 
ffie \\ = T - Thc ii- 

Lastly, to provide measurements that can be used to model 
the helium abundance in the envelope of the stars, following B14 
we extract the amplitude of the signal at Vmax from Eq. ([^, ob¬ 
taining 

Amax — AVjYiax OXp {—2b CL>Ynsix) ’ (4) 


with tUmax = ^TTVniax^ and v^ax derived from the background fit 
done by Cl5. Following [Verma et al.[ ( [2Q14] ), we also derive the 
amplitude of the signal in the frequencies, Ahc, given as 


Ahc = 


4cos2 [27rTHeII<Ay)] ’ 


(5) 


where (Ay) is the same mean large frequency separation used 
to calculate id and thc ii is the same acoustic depth used in 
Eq. 0. The adoption of Ahc to retrieve the helium content in 


2 
























































E. Corsaro et al.: High-precision acoustic helium signatures in 18 low-mass low-luminosity red giants 


the envelope is preferred since this value is not influenced by 


the location of the glitch (see e.g. 

Mazumdar & Antia[ [20011 

[Mazumdar et al.[[2014| [Verma et a] 

.[ 2014[). Eor clarity to the 

reader, we stress that Amax is derivec 

from a Bayesian approach 


by using the same sampling of the posterior probability distribu¬ 
tion obtained by Diamonds for the free para meters of the glitch 
model (see also [Corsaro & De Ridder|2()T4^ Fig. 7, for an anal¬ 
ogous case presenting the sampling from Diamonds for a com¬ 
bination of different inferred parameters). The parameters id 
and Ahc, simply follow from their definitions presented above, 
by using both the value (Av) computed from the radial mode 
frequencies presented by C15 for each star of our sample, and 
the estimated model parameters of Eq. 0 (see Sect. for more 
details). 


3. Results 


The results for Av^ (y), and (v) for the star 

KIC 12008916 are shown in Fig. (top and bottom panel, 
respectively), and can be found in Appendix for all the 
other LRGs, together with the tables with the individual mea¬ 
surements of the radial angular frequencies and of the 
corresponding second angular frequency differences (^) 

used in this work. We discarded KIC 10123207 from the fit 
because of the low number of available measurements (four 
second differences only, one less than the minimum required 
to fit the model given by Eq. [^. To help the reader visualize 
the presence of the oscillatory signal in the first differences of 
the angular degrees I = 0,2,3 and in the second differences 
of the angular degrees I = 2,3, we included low-degree (3-4) 
polynomial fits. The l-cr uncertainties on the first and second 
differences derived from a standard error propagation of the 
uncertainties of the individual mode frequencies, following 
Eq. Q and Eq. respectively, are overlaid in each plot, 
though they are not visible in most of the cases because they are 
smaller than the size of the symbols used for the measurements. 
The uncertainties on all the measurements are listed in the 
corresponding tables in Appendix [A| for completeness. 

We find that all the stars analyzed have clear acoustic glitches 
due to the signature of the Hen zone in ^ = 0 and 2 modes 
up to the second differences. We can see the presence of acous¬ 
tic glitches also in 7” = 3 modes for most of the stars thanks 
to at least four different frequency measurements available. As 
mentioned in Sect. we find that the measurements for modes 
having angular degrees 7 = 2,3 often deviate fro m t hose o f the 
radial oscillations (see e.g. Figs. A.6 |A.8 A.IO and A.13 ). As 
also indicated by C15 for the case of the mode linewidths of the 
7 = 2 modes, this different behavior relies on the presence of 
both mixed quadrupole modes and rotational split components. 
When using an individual Lorentzian profile to fit the frequency 
region containing the oscillation peak, either an 7 = 2 or 3 mode, 
as done by Cl5, the effects mentioned before can significantly 
change the measured frequency of the peak. A reliable treatment 
of the mixed modes and of the rotational split components for 
7 = 2, 3 modes is however difficult due to the high proximity of 
the individual peaks. 

The model fit to the acoustic glitch signatures of 
KIC 12008916 is shown in the bottom panel of Fig. for the 
case of A 2 y„,o (^) (red line), and similarly for the other stars in 
Appendix All the estimated parameters of Eq. ^ are pro¬ 
vided in Table [2 The inset in the bottom panel of Eig. pro¬ 
vides a closer view of one of the measurements to visualize the 
precision-level achieved in the fit. In particular, we find that the 
model given by Eq. ^ yields a remarkable fit quality for most of 



Fig. 1. Acoustic glitches for KIC 12008916. Top panel, the first 
difference, Ay(y), Eq. Q- Blue squares represent values com¬ 
puted from 7 = 0 modes, while green triangles from 7 = 2 
modes, with polynomial fits (dashed lines with same color as 
symbols) overlaid to visualize the oscillatory trend. The dashed 
vertical line marks the position of y^ax for a reference. Bottom 
panel the second difference A 2 y (y), Eq. Q, with same symbol 
description as for the top panel. The solid red line indicates the 
fit to the 7 = 0 measurements given by Eq. with the estimated 
parameters listed in Tableas derived by Diamonds. The hori¬ 
zontal dotted line marks the offset level cliln), useful to visual¬ 
ize the amplitude of the signature. The inset shows a zoom in of 
one of the measurements to visualize the precision-level of the 
fit. 


the stars, with average uncertainties of ~2% for thc id and ~ 8 % 
for b. Eollowing the analysis presented by [Corsaro et al.|p013| ), 
we have obtained the weighted Caussian rms of the residuals, 
(Tnns (fisted in Table as well). Eor its computation we have 
adopted the weights w/ = ct/ being the uncertainties on 
the second frequency difference coming from those reported in 
Table [AT] and |A.2] The quantity crrms provides additional in¬ 
formation to the reader because it allows to compare the qual¬ 
ity of the fits between different stars, and to relate the precision 
achieved on the individual fits to the given uncertainties of the 
data points. We note that for all the fits presented in this work, the 
values for (Trms are remarkably small, ranging from 10“^ down to 
10“^ yuHz in the best cases, thus reaching in many cases the same 
order of the precision-level obtained on the individual frequen¬ 
cies of the radial modes. Eor a reference to the reader, in Table [T] 
we also provide the values for the total acoustic radius T, with its 
l-cr standard deviation, and the values of y^ax obtained by C15. 
The stars KIC 8475025, KIC 9145955, KIC 10200377, and 


KIC 11913545, each show a component at high frequency that 
is not properly predicted by the adopted model. This mainly re¬ 
lies on some residual frequency dependence of the second dif¬ 
ferences that becomes more pronounced toward the wings of 
the region containing the oscillations. However, we note that 
the measurements at higher frequencies all have larger error- 
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bars (up to 10 times) with respect to the others, because of the 
larger linewidths of the peaks occurring at high-frequency (see 
C15 for more details). The fits derived, except for KIC 8475025, 
are therefore not significantly affected by the measurements of 
the second differences falling at high-frequency, whereas they 
are almost entirely constrained by those close to Vmax- This 
is inspected by refitting the glitch model without the highest- 
frequency measurement of the second difference (showing the 
deviating behavior), hence noticing that the new estimated free 
parameters of the model lay well within the uncertainties of 
those reported in Table For KIC 8475025 however, we find 
that the fit is unstable due to the large deviation (more than 
0.2/iHz) of the second dif feren ce measurement falling at the 
highest frequency (see Fig. |A.7| ). This is because the measured 
oscillation frequency of the highest-frequency radial mode is 
likely affected by additional sources such as mixed modes and 
rotational split components arising from the neighbor i = 2 
mode, which are enhanced by the large mode linewidths (see 
C15 for more details). To stabilize the fit for this star and pro¬ 
vide estimates of the model parameters that are comparable to 
the other stars in the sample, we have therefore chosen to dis¬ 
card the last measurement of the second difference for this par¬ 
ticular target. In the case of KIC 8366239, KIC 9267654, and 
KIC 10200377, the highest-frequency values are marked as not 
reliable (open symbols), according to the Bayesian peak signif¬ 
icance test done by C15, although they were included in the fit 
since they don’t produce any significant change in the results for 
the same reasons discussed above. 

The measurements of the acoustic radius ^He ii, the amplitude 
of the signal A^ax f^om Eq. and the corresponding charac¬ 
teristic width h, are shown in Fig. (top, middle and bottom 
panels, respectively) for all the stars of the sample. We note 
that while the model parameters (Eq. and their correspond¬ 
ing 68 % Bayesian credible intervals are estimated by means 
of Diamonds (see Table [T] and [Corsaro & De Ridd^|2014| for 
more details on the derivation of the Bayesian uncertainties), the 
uncertainties for the additional parameters ^He id ^max, and Ahc, 
were obtained in a subsequent step. In particular for A^ax, we 
have used the same sampling of the posterior probability distri¬ 
bution obtained by Diamonds, hence we have derived the me¬ 
dian and the corresponding 68 % Bayesian credible intervals di¬ 
rectly from the marginal probability distributions of Amax- For 
^He II and Ahc, the uncertainties follow from those of the acous¬ 
tic depth through the definition of the acoustic radius, and by a 
rescaling of the uncertainties on Amax through Eq. respec¬ 
tively. All the resulting values are listed in Table as well. We 
note that the precision obtained on our measurements of the 
acoustic radii o f the Hen zones is about 10 times higher than 
that obtained by jMiglio et al.| ( [20TQ| ) using CoRoT data. In addi¬ 
tion, all the values match those predicted by B14 along the entire 
range of Vmax investigated, showing a clear increasing trend to¬ 
ward lower Vmax, as expected for more evolved stages of the evo¬ 
lution in the RGB. The derived amplitudes in frequency, Ahc, 
are within the range 0.08-0.33 /iHz, and are varying from star to 
star with uncertainties on average around ~6 %, thus opening the 
possibility to study the He abundance by direct comparison with 
stellar models. 


4. Conclusions 

By exploiting the set of individual mode frequencies extracted 
by C15 for a sample of 19 LRGs with a precision level up to 
10“^ yuHz, we computed the first differences, Eq. ([^, and the sec¬ 
ond differences, Eq. ([^, for presenting the results on the acous- 



120 140 160 180 

I^max (mHz) 


Fig. 2. Acoustic radius of the Heii zone (top panel) and cor¬ 
responding amplitude of the oscillatory signal at Vmax (middle 
panel), and characteristic width b (bottom panel), as a function 
of Vmax for all the stars of the sample. The 68 % Bayesian uncer¬ 
tainties listed in Table [T] are overlaid for all the measurements. 


tic glitches of these stars. In this work, we have shown that the 
acoustic glitches are remarkably clear for all the red giants of the 
sample, and for both £ = 0 and £ = 2 modes, up to the second 
difference (where five to seven different measurements are avail¬ 
able for each star, except KIC 10123207 that instead has only 
four and was not considered in the analysis), with many cases 
involving £ = 3 modes as well. 

By adopting the model for the second differences introduced 
by [Houdek & Gough] ( |2007| ), Eq. ([^, we have extracted the 
acoustic depth, the characteristic width and amplitude of the sig¬ 
nal of all the Hen zones of the stars in our sample (see Table [T]) 
with an unprecedented level of detail for red giant stars (on av¬ 
erage ~2 % for the acoustic depths, ~8 % for the characteristic 
widths, and ~6 % for the amplitudes of the glitch signal in both 
second difference and frequency), improved by about one order 
of magnitude with respect to existing measurements of acous¬ 
tic depths from the literature. These values, reflecting the high- 
precision obtained on the individual frequency measurements of 
the radial modes, confirm the theoretical predictions done by 
B14 in the same range of Vmax- We also stress that the glitch 
model given by [Houdek & Goughj ( [2007 [ ) is able to predict the 
observations very exhaustively for most of the stars (well within 
the given uncertainties of the measurements for most of the data 
points available, as shown in Eig.l^and in the other figures in the 
Appendix). This is also supported by our computation of (Tnns, 
listed in Table which are on the same precision-level of that 
given by the measurements of the second difference for most of 
the stars analyzed. 

Einally, the set of values for A^ax and Ahc derived in this 
work, the latter not being influenced by the position of the glitch, 
coupled with the high precision achieved, will be useful to inves¬ 
tigate the helium content in the envelope of the stars, and possi- 
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Table 1. Median values with corresponding 68.3 % shortest credible intervals as derived by Diamonds for the free parameters of 
the model to fit the acoustic glitch signatures in the second differences, Eq. ([^, and for the acoustic radius of the He ii zone, the total 
acoustic radius of the star, the amplitude at v^ax in the second difference and in frequency (Eqs.|^an d0 respectively). 


KICID 
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A 

b 

<!> 

c 

^Hell 

T 

A 
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0.186:!;!11 
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1176:* 

4 1 +0.5 
^•^-0.4 

n 7^+0.08 
^•'•^-0.06 

27402:]| 

36634:110 

0 163+0-013 

u.iUC»_0.014 

0.083:!;*] 

182.04 

0.011 

10200377 

15660:g« 


1517:61 

0 9+0.5 
^•^-0.6 

0-76:6;* 

24327:™ 

39988+11 

0 147+0.017 
-0.024 

0 330+0-038 

142.52 

0.029 

10257278 

10565:^^1 

0.1106:«;«18‘ 

1639:g 

1 5+0-3 

0.83:6;05 

30551:]‘] 

41116:1* 

0 144+0.007 
-0.008 

0 075+0-004 

U.U/ J-0 004 

149.47 

0.007 

11353313 

16819!^* 

0.0104:«;««28 

1255:1“ 

1 4+0.3 
^•^-0.3 

0.71:6;04 

29749+1] 

46568:111 

0.178:6;!11 

0 240+0-020 
^•^^^-0.015 

126.46 

0.006 

11913545 

15218:228 


1412:81 

o 4+0.3 
'^•^-0.3 

0-70:6;62 

34222:*] 

49440*11 

0-164:!;*! 

0 127+0-005 
-0.005 

117.16 

0.005 

11968334 

13429:2«4 

-0.0015 

987+68 

2 6+0-3 

^•^-0.3 

0 72+0-04 
' -0.04 

30573+^6 

44003:111 

n 997+0.007 
vj.z.z.O-0 009 

0 169+0005 
^•^'^^-0.007 

141.43 

0.007 

12008916 

12492:286 

0 0113+^-^^^^ 

987:62 

2 0+0-3 

^•^-0.3 

n 77+0.04 
'-^-0.05 

26505:^84 

38998:1* 

0 247+0.007 

0 91 9+0.006 
u.z.iz,-ooio 

161.92 

0.013 


Notes. The parameters refer to the angular measurements of the second differences, A 2 du„,o (D» and the corresponding frequencies CL>n,o, of the 
radial modes only. 

Notes. The last two columns provide the reference values for Vmax (provided by C15) and the weighted Gaussian rms of the residuals, as described 
in Sect.|3] 


bly contribute to study scenarios of helium enrichment in low- 
mass stars (e.g. see |Gratton et al.|2012| and references therein). 
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Appendix A: Results for the fitting of the acoustic 
signatures 

All the individual angular frequency measurements and corre- 
spondi ng se con d diff erences of the radial modes are listed in the 
Tables |A.1| and |A.2| The results for the first differences Av^ (v) 
and second differ ences A^v^) of all the LRGs, are shown 

similarly to that provided for 


in the Figs, from [AT] to A. 17 
KIC 12008916inFTg:|2 
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with yellow star from i = 2 modes and corresponding polyno¬ 
mial fit with same color. 



Frequency (yuHz) 

Fig.A.2. Same description as in Fig. but for KIC 6117517, 
with yellow star from 1 = 2 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 



with yellow star from 1 = 2 modes and corresponding polyno¬ 
mial fit with same color. 
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Fig. A.5. Same description as in Fig.[^but for KIC 7619745. Fig. A.7. Same description as in Fig.[^but for KIC 8475025. 



Open symbols represent measurements that used modes with de¬ 
tection probability under the threshold suggested by Cl5. 



Frequency (yuHz) 


Fig.A.8. Same description as in Fig. but for KIC 8718745, 
with yellow star from i = 3 modes and corresponding polyno¬ 
mial fit with same color. 
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with yellow star from ^ = 3 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 



Fig. A. 10. Same description as in Fig. 0 but for KIC 9267654. 
Open symbols represent measurements that used modes with de¬ 
tection probability under the threshold suggested by Cl5. 



Fig.A.ll. Same description as in Fig. but for KIC 9475697, 
with yellow star from ^ = 3 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 



Fig. A.12. Same description as in Fig.[^but for KIC 9882316. 
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mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 



Frequency (yuHz) 

Fig. A.14. Same description as in Fig.[^but for KIC 10257278, 
with yellow star from 7” = 3 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 



with yellow star from 7 = 3 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 
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Fig. A.17. Same description as in Fig.[^but for KIC 11968334, 
with yellow star from ^ = 3 modes and corresponding polyno¬ 
mial fit with same color. Open symbols represent measurements 
that used modes with detection probability under the threshold 
suggested by C15. 
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Table A.l. Median values with corresponding 68.3 % shortest 
credible intervals for the radial angular frequencies and cor¬ 
responding second angular frequency differences A 2 tu„,o with 1- 
(T error bars, as described in Sect.|^ for all the stars of the sam¬ 
ple. 
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/ j 03 

1 '^n+O.Ol 
^•^'-0.07 
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Table A.2. Table A.l continued. 
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